Virological characteristics of the SARS-CoV-2 JN.1 variant (preprint)

The SARS-CoV-2 BA.2.86 lineage, first identified in August 2023, is phylogenetically distinct from the currently circulating SARS-CoV-2 Omicron XBB lineages, including EG.5.1 and HK.3. Comparing to XBB and BA.2, BA.2.86 carries more than 30 mutations in the spike (S) protein, indicating a high potential for immune evasion. BA.2.86 has evolved and its descendant, JN.1 (BA.2.86.1.1), emerged in late 2023. JN.1 harbors S:L455S and three mutations in non-S proteins. S:L455S is a hallmark mutation of JN.1: we have recently shown that HK.3 and other "FLip" variants carry S:L455F, which contributes to increased transmissibility and immune escape ability compared to the parental EG.5.1 variant. Here, we investigated the virological properties of JN.1.

Antiviral humoral immunity against SARS-CoV-2 Omicron subvariants induced by XBB.1.5 monovalent vaccine in infection-naive and XBB-infected individuals (preprint)

To control infection with SARS-CoV-2 Omicron XBB subvariants, the XBB.1.5 monovalent mRNA vaccine has been available since September 2023. However, we have found that natural infection with XBB subvariants, including XBB.1.5, does not efficiently induce humoral immunity against the infecting XBB subvariants. These observations raise the possibility that the XBB.1.5 monovalent vaccine may not be able to efficiently induce humoral immunity against emerging SARS-CoV-2 variants, including a variety of XBB subvariants (XBB.1.5, XBB.1.16, XBB.2.3, EG.5.1 and HK.3) as well as BA.2.86. To address this possibility, we collected two types of sera from individuals vaccinated with the XBB.1.5 vaccine; those who had not been previously infected with SARS-CoV-2 and those who had been infected with XBB subvariants prior to XBB.1.5 vaccination. We collected sera before and 3-4 weeks after vaccination, and then performed a neutralization assay using these sera and pseudoviruses.

Virological characteristics of the SARS-CoV-2 Omicron HK.3 variant harboring the "FLip" substitution (preprint)

In November 2023, SARS-CoV-2 XBB descendants, including EG.5.1 (XBB.1.9.2.5.1), the currently predominant lineage, are circulating worldwide according to Nextstrain. EG.5.1 has a characteristic amino acid substitution in the spike protein (S), S:F456L, which contributes to its escape from humoral immunity. EG.5.1 has further evolved, and its descendant lineage harboring S:L455F (i.e., EG.5.1+S:L455F) emerged and was named HK.3 (XBB.1.9.2.5.1.1.3). HK.3 was initially discovered in East Asia and is rapidly spreading worldwide. Notably, the XBB subvariants bearing both S:L455F and S:F456L substitutions, including HK.3, are called the "FLip" variants. These FLip variants, such as JG.3 (XBB.1.9.2.5.1.3.3), JF.1 (XBB.1.16.6.1) and GK.3 (XBB.1.5.70.3), have emerged convergently, suggesting that the acquisition of these two substitutions confers a growth advantage to XBB in the human population. Here, we investigated the virological properties of HK.3 as a representative of the FLip variants.

Virological characteristics of the SARS-CoV-2 BA.2.86 variant (preprint)

In late 2023, a lineage of SARS-CoV-2 emerged and was named the BA.2.86 variant. BA.2.86 is phylogenetically distinct from other Omicron sublineages identified so far, displaying an accumulation of over 30 amino acid mutations in its spike protein. Here, we performed multiscale investigations to reveal the virological characteristics of the BA.2.86 variant. Our epidemic dynamics modeling suggested that the relative reproduction number of BA.2.86 is significantly higher than that of EG.5.1. Experimental studies showed that four clinically-available antivirals were effective against BA.2.86. Although the fusogenicity of BA.2.86 spike is similar to that of the parental BA.2 spike, the intrinsic pathogenicity of BA.2.86 in hamsters was significantly lower than that of BA.2. Since the growth kinetics of BA.2.86 is significantly lower than that of BA.2 in both in vitro cell cultures and in vivo, it is suggested that the attenuated pathogenicity of BA.2.86 is due to its decreased replication capacity.

Virological characteristics of the SARS-CoV-2 Omicron EG.5.1 variant (preprint)

In middle-late 2023, a sublineage of SARS-CoV-2 Omicron XBB, EG.5.1 (a progeny of XBB.1.9.2), is spreading rapidly around the world. Here, we performed multiscale investigations to reveal virological features of newly emerging EG.5.1 variant. Our phylogenetic-epidemic dynamics modeling suggested that two hallmark substitutions of EG.5.1, S:F456L and ORF9b:I5T, are critical to the increased viral fitness. Experimental investigations addressing the growth kinetics, sensitivity to clinically available antivirals, fusogenicity and pathogenicity of EG.5.1 suggested that the virological features of EG.5.1 is comparable to that of XBB.1.5. However, the cryo-electron microscopy reveals the structural difference between the spike proteins of EG.5.1 and XBB.1.5. We further assessed the impact of ORF9b:I5T on viral features, but it was almost negligible at least in our experimental setup. Our multiscale investigations provide the knowledge for understanding of the evolution trait of newly emerging pathogenic viruses in the human population.

Transmissibility, infectivity, and immune resistance of the SARS-CoV-2 BA.2.86 variant (preprint)

In September 2023, the SARS-CoV-2 XBB descendants, such as XBB.1.5 and EG.5.1 (originally XBB.1.9.2.5.1), are predominantly circulating worldwide. Unexpectedly, however, a lineage distinct from XBB was identified and named BA.2.86 on August 14, 2023. Notably, BA.2.86 bears more than 30 mutations in the spike (S) protein when compared to XBB and the parental BA.2, and many of them are assumed to be associated with immune evasion. Although the number of reported cases is low (68 sequences have been reported as of 7 September 2023), BA.2.86 has been detected in several continents (Europe, North America and Africa), suggesting that this variant may be spreading silently worldwide. On 17 August 2023, the WHO designated BA.2.86 as a variant under monitoring. Here we show evidence suggesting that BA.2.86 potentially has greater fitness than current circulating XBB variants including EG.5.1. The pseudovirus assay showed that the infectivity of BA.2.86 was significantly lower than that of B.1.1 and EG.5.1, suggesting that the increased fitness of BA.2.86 is not due to the increased infectivity. We then performed a neutralization assay using XBB breakthrough infection sera to address whether BA.2.86 evades the antiviral effect of the humoral immunity induced XBB subvariants. The 50% neutralization titer of XBB BTI sera against BA.2.86 was significantly (1.4-fold) lower than those against EG.5.1. The sera obtained from individuals vaccinated with 3rd-dose monovalent, 4th-dose monovalent, 4th-dose BA.1 bivalent, and 4th-dose BA.5 bivalent mRNA vaccines exhibited very little or no antiviral effects against BA.2.86. Moreover, the three monoclonal antibodies (Bebtelovimab, Sotrovimab and Tixagevimab), which worked against the parental BA.2, did not exhibit antiviral effects against BA.2.86. These results suggest that BA.2.86 is one of the most highly immune evasive variants ever.

Virological characteristics of the SARS-CoV-2 XBB.1.5 variant (preprint)

Circulation of SARS-CoV-2 Omicron XBB has resulted in the emergence of XBB.1.5, a new Variant of Interest. Our phylogenetic analysis suggests that XBB.1.5 evolved from XBB.1 by acquiring the F486P spike (S) mutation, subsequent to the acquisition of a nonsense mutation in ORF8. Neutralization assays showed similar abilities of immune escape between XBB.1.5 and XBB.1. We determined the structural basis for the interaction between human ACE2 and the S protein of XBB.1.5, showing similar overall structures between the S proteins of XBB.1 and XBB.1.5. The intrinsic pathogenicity of XBB.1.5 in hamsters is lower than that of XBB.1. Importantly, we found that the ORF8 nonsense mutation of XBB.1.5 resulted in impairment of MHC expression. In vivo experiments using recombinant viruses revealed that the XBB.1.5 mutations are involved with reduced virulence of XBB.1.5. Together, these data suggest that the mutations in ORF8 and S could enhance spreading of XBB.1.5 in humans.

Antiviral efficacy of the SARS-CoV-2 XBB breakthrough infection sera against Omicron subvariants including EG.5 (preprint)

As of July 2023, EG.5.1 (a.k.a. XBB.1.9.2.5.1), a XBB subvariant bearing the S:Q52H and S:F456L substitutions, alongside the S:F486P substitution (Figure S1A), has rapidly spread in some countries. On July 19, 2023, the WHO classified EG.5 as a variant under monitoring. First, we showed that EG.5.1 exhibits a higher effective reproduction number compared with XBB.1.5, XBB.1.16, and its parental lineage (XBB.1.9.2), suggesting that EG.5.1 will spread globally and outcompete these XBB subvariants in the near future. We then addressed whether EG.5.1 evades from the antiviral effect of the humoral immunity induced by breakthrough infection (BTI) of XBB subvariants and performed a neutralization assay using XBB BTI sera. However, the 50% neutralization titer (NT50) of XBB BTI sera against EG.5.1 was comparable to those against XBB.1.5/1.9.2 and XBB.1.16. Moreover, the sensitivity of EG.5.1 to convalescent sera of XBB.1- and XBB.1.5-infected hamsters was similar to those of XBB.1.5/1.9 and XBB.1.16. These results suggest that the increased Re of EG.5.1 is attributed to neither increased infectivity nor immune evasion from XBB BTI, and the emergence and spread of EG.5 is driven by the other pressures. We previously demonstrated that Omicron BTI cannot efficiently induce antiviral humoral immunity against the variant infected. In fact, the NT50s of the BTI sera of Omicron BA.1, BA.2, and BA.5 against the variant infected were 3.0-, 2.2-, and 3.4-fold lower than that against the ancestral B.1.1 variant, respectively. However, strikingly, we found that the NT50 of the BTI sera of XBB1.5/1.9 and XBB.1.16 against the variant infected were 8.7- and 8.3-fold lower than that against the B.1.1 variant. These results suggest that XBB BTI cannot efficiently induce antiviral humoral immunity against XBB subvariants.

Determination of the factors responsible for host tropism of SARS-CoV-2-related bat coronaviruses (preprint)

Differences in host ACE2 genes may affect the host range of SARS-CoV-2-related coronaviruses (SC2r-CoVs) and further determine the tropism of host ACE2 for the infection receptor. However, the factor(s) responsible for determining the host tropism of SC2r-CoVs, which may in part be determined by the tropism of host ACE2 usage, remains unclear. Here, we use the pseudoviruses with the spike proteins of two Laotian SC2r-CoVs, BANAL-20-236 and BANAL-20-52, and the cells expressing ACE2 proteins of eight different Rhinolophus bat species, and show that these two spikes have different tropisms for Rhinolophus bat ACE2. Through structural analysis and cell culture experiments, we demonstrate that this tropism is determined by residue 493 of the spike and residues 31 and 35 of ACE2. Our results suggest that SC2r-CoVs exhibit differential ACE2 tropism, which may be driven by adaptation to different Rhinolophus bat ACE2 proteins.

Virological characteristics of the SARS-CoV-2 Omicron XBB.1.16 variant (preprint)

At the end of March 2023, XBB.1.16, a SARS-CoV-2 omicron XBB subvariant, emerged and was detected in various countries. Compared to XBB.1.5, XBB.1.16 has two substitutions in the S protein: E180V is in the N-terminal domain, and T478K in the receptor-binding domain (RBD). We first show that XBB.1.16 had an effective reproductive number (Re) that was 1.27- and 1.17-fold higher than the parental XBB.1 and XBB.1.5, respectively, suggesting that XBB.1.16 will spread worldwide in the near future. In fact, the WHO classified XBB.1.16 as a variant under monitoring on March 30, 2023. Neutralization assays demonstrated the robust resistance of XBB.1.16 to breakthrough infection sera of BA.2 (18-fold versus B.1.1) and BA.5 (37-fold versus B.1.1). We then used six clinically-available monoclonal antibodies and showed that only sotrovimab exhibits antiviral activity against XBB subvariants, including XBB.1.16. Our results suggest that, similar to XBB.1 and XBB.1.5, XBB.1.16 is robustly resistant to a variety of anti-SARS-CoV-2 antibodies. Our multiscale investigations suggest that XBB.1.16 that XBB.1.16 has a greater growth advantage in the human population compared to XBB.1 and XBB.1.5, while the ability of XBB.1.16 to exhibit profound immune evasion is comparable to XBB.1 and XBB.1.5. The increased fitness of XBB.1.16 may be due to (1) different antigenicity than XBB.1.5; and/or (2) the mutations in the non-S viral protein(s) that may contribute to increased viral growth efficiency.

SARS-CoV-2 ORF3c suppresses immune activation by inhibiting innate sensing (preprint)

SARS-CoV-2 proteins are translated from subgenomic RNAs (sgRNAs). While most of these sgRNAs are monocistronic, some viral mRNAs encode more than one protein. For example, the ORF3a sgRNA also encodes ORF3c, an enigmatic 41-amino acid peptide. Here, we show that ORF3c suppresses RIG-I- and MDA5-mediated immune activation and interacts with the signaling adaptor MAVS. In line with this, ORF3c inhibits IFN-{beta} induction. This immunosuppressive activity of ORF3c is conserved among members of the subgenus sarbecovirus, including SARS-CoV and coronaviruses isolated from bats. Notably, however, the SARS-CoV-2 delta and kappa variants harbor premature stop codons in ORF3c demonstrating that this reading frame is not essential for efficient viral replication in vivo. In agreement with this, disruption of ORF3c did not significantly affect SARS-CoV-2 replication in CaCo-2 or CaLu-3 cells. In summary, we here identify ORF3c as an immune evasion factor that suppresses IFN-{beta} induction, but is dispensable for efficient replication of SARS-CoV-2.

Enhanced transmissibility, infectivity and immune resistance of the SARS-CoV-2 Omicron XBB.1.5 variant (preprint)

In 2022, we have elucidated the characteristics of a variety of newly emerging SARS-CoV-2 Omicron subvariants. At the end of 2022, the XBB.1.5 variant, an descendant of XBB.1 that acquired the S:F486P substitution, emerged and is rapidly spreading in the USA and is the latest variant of concern. Although the features of XBB.1.5 was already reported by another group as a preprint, we think multiple and independent evaluations important, and these reports are crucial for sustained global health. In this study, our epidemic dynamics analysis revealed that the relative effective reproduction number (Re) of XBB.1.5 is more than 1.2-fold greater than that of the parental XBB.1, and XBB.1.5 is outcompeting BQ.1.1, the predominant lineage in the USA as of December 2022. Our data suggest that XBB.1.5 will rapidly spread worldwide in the near future. Yeast surface display assay and pseudovirus assay respectively showed that the ACE2 binding affinity and infectivity of XBB.1.5 is 4.3-fold and 3.3-fold higher than those of XBB.1, respectively. Moreover, neutralization assay revealed that XBB.1.5 is robustly resistant to BA.2 breakthrough infection sera (41-fold versus B.1.1, 20-fold versus BA.2) and BA.5 breakthrough infection sera (32-fold versus B.1.1, 9.5-fold versus BA.5), respectively. Because the immune resistance of XBB.1.5 is comparable to that of XBB.1, our results suggest that XBB.1.5 is the most successful XBB lineage as of January 2023 by acquiring the S:F486P substitution to augment ACE2 binding affinity without losing remarkable immune resistance, which leads to greater transmissibility.

Virological characteristics of the SARS-CoV-2 XBB variant derived from recombination of two Omicron subvariants (preprint)

In late 2022, the SARS-CoV-2 Omicron subvariants have highly diversified, and XBB is spreading rapidly around the world. Our phylogenetic analyses suggested that XBB emerged by recombination of two co-circulating BA.2 lineages, BJ.1 and BM.1.1.1 (a progeny of BA.2.75), during the summer of 2022 around India. In vitro experiments revealed that XBB is the most profoundly resistant variant to BA.2/5 breakthrough infection sera ever and is more fusogenic than BA.2.75. Notably, the recombination breakpoint is located in the receptor-binding domain of spike, and each region of recombined spike conferred immune evasion and augmented fusogenicity to the XBB spike. Finally, the intrinsic pathogenicity of XBB in hamsters is comparable to or even lower than that of BA.2.75. Our multiscale investigation provided evidence suggesting that XBB is the first documented SARS-CoV-2 variant increasing its fitness through recombination rather than single mutations.

Convergent evolution of the SARS-CoV-2 Omicron subvariants leading to the emergence of BQ.1.1 variant (preprint)

In late 2022, although the SARS-CoV-2 Omicron subvariants have highly diversified, some lineages have convergently acquired amino acid substitutions at five critical residues in the spike protein. Here, we illuminated the evolutionary rules underlying the convergent evolution of Omicron subvariants and the properties of one of the latest lineages of concern, BQ.1.1. Our phylogenetic and epidemic dynamics analyses suggest that Omicron subvariants independently increased their viral fitness by acquiring the convergent substitutions. Particularly, BQ.1.1, which harbors all five convergent substitutions, shows the highest fitness among the viruses investigated. Neutralization assays show that BQ.1.1 is more resistant to breakthrough BA.2/5 infection sera than BA.5. The BQ.1.1 spike exhibits enhanced binding affinity to human ACE2 receptor and greater fusogenicity than the BA.5 spike. However, the pathogenicity of BQ.1.1 in hamsters is comparable to or even lower than that of BA.5. Our multiscale investigations provide insights into the evolutionary trajectory of Omicron subvariants.

Virological characteristics of the SARS-CoV-2 Omicron BA.2.75 (preprint)

SARS-CoV-2 Omicron BA.2.75 emerged in May 2022. BA.2.75 is a BA.2 descendant but is phylogenetically different from BA.5, the currently predominant BA.2 descendant. Here, we showed that the effective reproduction number of BA.2.75 is greater than that of BA.5. While the sensitivity of BA.2.75 to vaccination- and BA.1/2 breakthrough infection-induced humoral immunity was comparable to that of BA.2, the immunogenicity of BA.2.75 was different from that of BA.2 and BA.5. Three clinically-available antiviral drugs were effective against BA.2.75. BA.2.75 spike exhibited a profound higher affinity to human ACE2 than BA.2 and BA.5 spikes. The fusogenicity, growth efficiency in human alveolar epithelial cells, and intrinsic pathogenicity in hamsters of BA.2.75 were comparable to those of BA.5 but were greater than those of BA.2. Our multiscale investigations suggest that BA.2.75 acquired virological properties independently of BA.5, and the potential risk of BA.2.75 to global health is greater than that of BA.5.

Comparative pathogenicity of SARS-CoV-2 Omicron subvariants including BA.1, BA.2, and BA.5 (preprint)

Unremitting emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants imposes us to continuous control measurement. Given the rapid spread, new Omicron subvariant named BA.5 is urgently required for characterization. Here we analyzed BA.5 with the other Omicron variants BA.1, BA.2, and ancestral B.1.1 comprehensively. Although in vitro growth kinetics of BA.5 is comparable among the Omicron subvariants, BA.5 become much more fusogenic than BA.1 and BA.2. The airway-on-a-chip analysis showed that the ability of BA.5 to disrupt the respiratory epithelial and endothelial barriers is enhanced among Omicron subvariants. Furthermore, in our hamster model, in vivo replication of BA.5 is comparable with that of the other Omicrons and less than that of the ancestral B.1.1. Importantly, inflammatory response against BA.5 is strong compared with BA.1 and BA.2. Our data suggest that BA.5 is still low pathogenic compared to ancestral strain but evolved to induce enhanced inflammation when compared to prior Omicron subvariants.

Neutralization sensitivity of Omicron BA.2.75 to therapeutic monoclonal antibodies (preprint)

Since the end of 2021, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant outcompeted other variants and took over the world. After the emergence of original Omicron BA.1, Omicron BA.2 subvariant emerged and outcompeted BA.1. As of July 2022, some BA.2 subvariants, including BA.2.12.1, BA.4 and BA.5, emerged in multiple countries and begun outcompeting original BA.2. Moreover, a novel BA.2 subvariant, BA.2.75, was detected in eight countries including India at the end of June 2022, and preliminary investigations suggest that BA.2.75 is more transmissible over the other BA.2 subvariants. On July 7, 2022, the WHO classified BA.2.75 as a variant-of-concern lineage under monitoring. We have recently demonstrated that BA.4/5 is highly resistant to a therapeutic monoclonal antibody, cilgavimab, than BA.2. The resistance of SARS-CoV-2 variants to therapeutic antibodies can be attributed to the mutations in the viral spike protein. Compared to the BA.2 spike, BA.2.12.1 and BA.4/5 respectively bear two and four mutations in their spike proteins. On the other hand, the majority of BA.2.75 spike bears nine substitutions. The fact that the mutation number in the BA.2.75 spike is larger than those in the BA.4/5 spike raises the possibility that the BA.2.75 spike significantly reduces sensitivity towards therapeutic monoclonal antibodies than BA.2 and BA.4/5. In this study, we generated pseudoviruses harboring the spike proteins of BA.2.75, BA.4/5 and BA.2 and evaluated the efficacy of ten therapeutic monoclonal antibodies and three antibody cocktails against BA.2.75.

Virological characteristics of the novel SARS-CoV-2 Omicron variants including BA.2.12.1, BA.4 and BA.5 (preprint)

After the global spread of SARS-CoV-2 Omicron BA.2 lineage, some BA.2-related variants that acquire mutations in the L452 residue of spike protein, such as BA.2.9.1 and BA.2.13 (L452M), BA.2.12.1 (L452Q), and BA.2.11, BA.4 and BA.5 (L452R), emerged in multiple countries. Our statistical analysis showed that the effective reproduction numbers of these L452R/M/Q-bearing BA.2-related Omicron variants are greater than that of the original BA.2. Neutralization experiments revealed that the immunity induced by BA.1 and BA.2 infections is less effective against BA.4/5. Cell culture experiments showed that BA.2.12.1 and BA.4/5 replicate more efficiently in human alveolar epithelial cells than BA.2, and particularly, BA.4/5 is more fusogenic than BA.2. Furthermore, infection experiments using hamsters indicated that BA.4/5 is more pathogenic than BA.2. Altogether, our multiscale investigations suggest that the risk of L452R/M/Q-bearing BA.2-related Omicron variants, particularly BA.4 and BA.5, to global health is potentially greater than that of original BA.2.

Sensitivity of novel SARS-CoV-2 Omicron subvariants, BA.2.11, BA.2.12.1, BA.4 and BA.5 to therapeutic monoclonal antibodies (preprint)

As of May 2022, Omicron BA.2 variant is the most dominant variant in the world. Thereafter, Omicron subvariants have emerged and some of them began outcompeting BA.2 in multiple countries. For instance, Omicron BA.2.11, BA.2.12.1 and BA.4/5 subvariants are becoming dominant in France, the USA and South Africa, respectively. In this study, we evaluated the sensitivity of these new Omicron subvariants (BA.2.11, BA.2.12.1 and BA.4/5) to eight therapeutic monoclonal antibodies (bamlanivimab, bebtelovimab, casirivimab, cilgavimab, etesevimab, imdevimab, sotrovimab and tixagevimab). Notably, we showed that although cilgavimab is antiviral against BA.2, BA.4/5 exhibits higher resistance to this antibody compared to BA.2. Since mutations are accumulated in the spike proteins of newly emerging SARS-CoV-2 variants, we suggest the importance of rapid evaluation of the efficiency of therapeutic monoclonal antibodies against novel SARS-CoV-2 variants.

SARS-CoV-2 spike S375F mutation characterizes the Omicron BA.1 variant (preprint)

Recent studies have revealed the unique virological characteristics of Omicron, the newest SARS-CoV-2 variant of concern, such as pronounced resistance to vaccine-induced neutralizing antibodies, less efficient cleavage of the spike protein, and poor fusogenicity. However, it remains unclear which mutation(s) in the spike protein determine the virological characteristics of Omicron. Here, we show that the representative characteristics of the Omicron spike are determined by its receptor-binding domain. Interestingly, the molecular phylogenetic analysis revealed that the acquisition of the spike S375F mutation was closely associated with the explosive spread of Omicron in the human population. We further elucidate that the F375 residue forms an interprotomer pi-pi interaction with the H505 residue in another protomer in the spike trimer, which confers the attenuated spike cleavage efficiency and fusogenicity of Omicron. Our data shed light on the evolutionary events underlying Omicron emergence at the molecular level.